Capacities of Quantum Channels for Massive Bosons and Fermions 
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We consider the capacity of classical information transfer for noiseless quantum channels carrying 
a finite average number of massive bosons and fermions. The maximum capacity is attained by 
transferring the Fock states generated from the grand-canonical ensemble. Interestingly, the chan- 
nel capacity for a Bose gas indicates the onset of a Bose-Einstein condensation, by changing its 
qualitative behavior at the criticality, while for a channel carrying weakly attractive fermions, it ex- 
hibits the signatures of Bardeen-Cooper-Schrieffer transition. We also show that for noninteracting 
particles, fermions are better carriers of information than bosons. 



A communication channel carrying classical informa- 
tion by using quantum states as the carriers of informa- 
tion, has been a subject of intensive studies. The fun- 
damental result in this respect is the "Holevo bound" 
(see also e.g. obtained more than 30 years 

ago, which gives the capacity of such channels. An es- 
sential message carried by the Holevo bound is that at 
most n bits (binary digits) of classical information can 
be sent via a quantum system of n distinguishable qubits 
(two-dimensional quantum systems). However in realistic 
channels, where the quantum system is usually of infinite 
dimensions, the Holevo bound predicts infinite capacities. 
In realistic channels, it is therefore important to give a 
physical constraint on the carriers of the information. 

Information carried over long distances usually em- 
ploys electromagnetic signals as carriers of information. 
Capacities of such channels have been studied quite ex- 
tensively (see e.g. 0, 0, @). In this case, the physi- 
cal constraint that is used to avoid the infinite capacity 
problem, is an energy constraint. Due to the form of 
the Holevo bound, the ensemble that maximizes the ca- 
pacity, turns out to be the canonical ensemble (or the 
microcanonical ensemble, depending on the type of the 
energy constraint) of statistical mechanics |5j- 

In recent experiments, it has been possible toproduce 
atomic waveguides in optical microstructures [7(, or on 
an atom chip that may serve as quantum channels 
of macroscopic (or at least mesoscopic) length scales. 
Channels carrying massive particles have possibly fasci- 
nating applications in quantum information processing. 
It is thus important to obtain their capacities. Since we 
are dealing now with massive information carriers, it is 
not enough to put an energy constraint only. Rather, it 
is natural to give a particle number constraint as well 
as an energy constraint. This, of course, hints at the 
grand-canonical ensemble (GCE) of statistical mechan- 
ics. Indeed, we show in this Letter, the ensemble that 
maximizes the capacity of noiseless channels that carry 
massive bosons or fermions, under particle number and 
energy constraints, is GCE. Note that massless photons, 
despite being bosons, do not exhibit Bose-Einstein con- 



densation (BEC), due to the lack of constraint on their 
number. Massive bosons, however, do exhibit BEC; we 
show in this Letter that the channel capacity of mas- 
sive bosons indicates the onset of BEC, by changing its 
behavior from being concave with respect to tempera- 
ture, to being convex. The bosons that we consider in 
this Letter are noninteracting. Noninteracting fermions, 
however, do not exhibit any phase transition. Interact- 
ing fermions, on the other hand, exhibit the Bardeen- 
Cooper-Schrieffer (BCS) transition, and as we show in 
this Letter, the capacity of interacting fermionic chan- 
nels exhibits the onset of such transition. We obtain our 
results by simulating a finite number of particles in the 
channel, and not the thermodynamical limit of an infinite 
number of particles. We also show that for a wide range 
of power law potentials, including the harmonic trap and 
the rectangular box, and for moderate and high tempera- 
tures, the fermions are better carriers of information than 
bosons, for the case of noninteracting particles. 

Suppose therefore that a sender (Alice) encodes the 
classical message i (occuring with probability pi) in the 
state Qi, and sends it to a receiver (Bob). The channel 
is noiseless, while Qi can be mixed. To obtain informa- 
tion about i, Bob performs a measurement M (on the 
ensemble £ = {pt, Qi}) to obtain the post-measurement 
ensemble {pi\ m , Qi\ m }i with probability q m . The infor- 
mation gained by this measurement can be quantified 
by the mutual information Im(£) between the index i 
and the measurement results m: Im(£) = H({Pi}) ~ 

E m ?™ ff ({ft|m})- Here H (in}) = -Ei r * lo S2 r i is the 
Shannon entropy of a probability distribution {r^}. The 

accessible information I acc (£) = maxM Zm" (£) is ob- 
tained by maximizing over all possible M. 

The Holevo bound gives a very useful upper bound 
on the accessible information for an arbitrary ensemble: 

Iacc(£) < x(£) = S(q) - J2iP* s (Qi)- Here Q = T,iPiQi, 
and S(n) = — tr (77 log 2 rj) is the von Neumann entropy 
of 77. In a noiseless environment, the capacity of such 
an information transfer is the maximum, over all input 
ensembles satisfying a given physical constraint, of the ac- 
cessible information. It is important to impose a physical 
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constraint on the input ensembles, as arbitrary encoding 
and decoding schemes are included in the Holevo bound. 
This has the consequence that the bound explodes for 
infinite dimensional systems: For an ensemble of pure 
states with average ensemble state g^, \ — S(g^,), which 
can be as large as log 2 d, where d is the dimension of 
the Hilbert space to which the ensemble belongs. If the 
pure states are orthogonal, I acc = \ = log 2 d, so that the 
capacity diverges along with its bound (see e.g. Q). 

To avoid this infinite capacity, one usually uses an en- 
ergy constraint for channels that carry photons (see e.g. 
0TIE01)- Suppose that the system is described by the 
Hamiltonian Ji. Then the average energy constraint on 
a communication channel that is sending the ensemble 
£ = {pi,Qi\, is tr(g>7Y) = E. Here g is the average en- 
semble state, and E is the average energy available to 
the system. The capacity Ce of such a channel is then 
the maximum of I a cc{£), over all ensembles, under the 
average energy constraint. Now I acc (£) < x(£) < S(g), 
and S(g) is maximized, under the same constraint, by the 
canonical ensemble (CE) corresponding to the Hamilto- 
nian H, and energy E (see e.g. :9|). Moreover, this is 
an ensemble of orthogonal pure states (the Fock states, 
or in other words number, states), so that Ce is also 
reached for this ensemble [3- The channel capacity de- 
pends solely on the average ensemble state, which in this 
optimal case is the canonical equilibrium (thermal) state 
g e i = exp(-/3H)/Z, where = l/k B T, with k B being 
the Boltzmann constant, and T the absolute tempera- 
ture. Z = tr(exp(— PTC)) is the partition function. For 
given E, T is given by E = — ^(log e Z). As a result 
of particle number nonconservation, noninteracting pho- 
tons and hence their capacity do not exhibit signatures 
of a condensation. However, effectively interacting pho- 
ton fluids and photon condensation effects are possible, 
in principle, by using nonlinear cavities, in which case 
the photons may acquire an effective mass (see e.g. [loj'). 

In the case of channels that carry massive particles, it 
is natural to impose the additional constraint of average 
particle number. Suppose that Alice prepares A particles 
in a trap, and transfers them to Bob. Let the trap have 
energy levels e^, and let ni be the average occupation 
number of the «-th level. Then the conservation of the 
average particle number reads = N, and the con- 

straint of a fixed average energy, for a given energy E, is 
X^j n i £ i = The channel capacity Ce,n of such a chan- 
nel is the maximum of I acc (£), over all ensembles that 
satisfy these two constraints. Under these constraints, 
the von Neumann entropy of the average state of the 
system is maximized by GCE. Again the ensemble ele- 
ments are pure and orthogonal (Fock states), whence the 
channel capacity is reached by the same ensemble. 

It is important to stress here that the channel capaci- 
ties that we derive in this Letter are all for the case of a 
given finite average number of particles in the trap, and 
not in the thermodynamic limit. This is because in a real 
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FIG. 1: The channel capacity, for the case of noninteracting 
bosons, plotted against T/T c . The lower curve is for the case 
of 100 bosons, in a 3D box with pbc. For 87 Rb atoms in such 
a box of volume 1/zm 3 , thermodynamical calculations predict 
T c ~ OAfiK. The upper curve is for 500 bosons in the same 
trap. The thermodynamical estimations of T c are higher than 
the corresponding values that we obtain (as indicated by the 
arrows) . 



implementation of such channels, this number is usually 
only at most moderately high. 

Noninteracting bosons. Here, n\ — 1/ (eP( £i ~ ti ) — 1), 
and the channel capacity (in bits) is given by (see e.g. 
E3) C b l N = -EiKlog2^-(l + ^)log 2 (l + n^)). 
Here fi b is the chemical potential. For given average par- 
ticle number A and absolute temperature T, one uses 
the energy constraint to find n b for that case. The en- 
ergy is then given by the average particle number con- 
straint, and the capacity by C b £ N . Let us consider the 
case when the trap is a 3D-box of volume L 3 with peri- 
odic boundary condition (pbc), so that the energy levels 

are lkl?-( n l + n v + n X7 n y ,n z = 0, ±1, . . .. Here 

m is the mass of the individual particles in the trap. 
In Fig. [J we plot C b £ N vs. T/T c , for different A. Here 

T c = (2vr ti 2 1 mk B ) (N / '2 .612 L 3 ) 2 ' ' 3 is the critical tempera- 
ture, as obtained in the thermodynamical (large N) limit. 
The capacity changes its shape from being concave to be- 
ing convex with respect to temperature, at the onset of 
a BEC. The thermodynamical estimation of T c is higher 
than our values of T c for different N, with the gap reduc- 
ing for growing N. Such indication of a gap has also been 
obtained previously (see e.g. [ll|). We have checked that 
the predicted approximate gap in Ref. [Tlj . is in agree- 
ment with our calculations for N = 1000. For lower N 
however, the prediction is no longer valid, as expected in 
Ref. [ijj- Note that the capacity increases with increas- 
ing A, and it has the same qualitative behavior for a 3D 
box without pbc, as well as for a harmonic trap. 

Fermions are better carriers of information than 
bosons. For spin-s noninteracting fermions, n\ — 
5 /( e /3(e,-M / ) + i) j w h ere g = 2s + 1, and {i f is 
the fermion chemical potential. The channel capac- 
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ity (in bits) is given by (see e.g. 12]) C E d N = 

~9J2i {jf lo S2 ^" + (1 - ^) log 2 (l - ^-)J . Again, for 
given N and T, one obtains pJ from average particle 
number conservation, which then gives the capacity. 

The bosons that we have considered in this Letter 
are spinless. To make a fair comparison of the capaci- 
ties, we consider "spinless" , i.e. polarized fermions with 
g=l. Let us start with bosons, and perform the high 
temperature expansion. First, we expand the fugacity 
z b = c^ in powers of N/S\, where Sk = J2 t e ~ kl3£i i 
and find the coefficients from the average particle num- 
ber conservation. We use this expansion to find an ex- 
pansion of the n b, s, which in turn is substituted in the 
formula for C b £ N . The same calculation is done for 
fermions. We perform the calculation up to the third 
order, and find that C b ^ N = (£f =1 a\ (N/SiY) lo S2 e + 
/3j(JV/Si)log 2 (JV/Si) + 'f3 b 2 (N/S 1 ) 2 \og 2 (N/S 1 ), whereas 
C f E d N = (E< 4 (N/S.T) log 2 e + (3{(N/S 1 )log 2 (N/S 1 ) + 
p((N/Siflog 2 (N/Si). The coefficients of first order 
perturbation are equal: a\ = a{ — S\ + D\. In 
the next order, they differ by a sign: a 2 — — a 2 — 
S 2 /2 - S 2 Di/Si + D 2 . The third order perturbation 
coefficients are again equal: a\ — — —3S 2 + S3/3 + 
25|/ 'Si+(2Si-SiS a )D 1 / Sf-2S 2 D 2 / Si+£ 3 . Also, fi\ = 
0l = -St, (3 b = 0, 0* = 2S 2 ; here D k = (3e t e~ k ^ . 

To find out the potentials and dimensions for which 
this perturbation technique is systematic, we consider 
uniform power law potentials, such as V — r 7 , in 
a d-dimensional Cartesian space of (x\, . . . , X4), with 
r = \Jx\ + . . . + x 2 d . We calculate Sk and D^, replac- 
ing sums by integrations with density of states p(e) = 
J d d pd d xS(p 2 /2m + V — e), the latter integration being 
over the phase space. One may then check that the tech- 
nique is systematic when I/7 + I/2 > 1/d. This includes, 
e.g., the 3D and 2D harmonic potentials, for which we 
also performed the summations directly, and obtained 
the same results. Note that = (3 2 < (3 2 , whereas 
a 2 = "5*2 (1 — P {e)a + (3 {s) 2 p) (where (s)a stands for the 
average energy with Boltzmann probabilities e~@ £i /Si) 
can be explicitly evaluated by using the density of states 
p(e). As a result, we obtain a 2 < < a 2 , implying that 
N < f° r a l ar g e range of sufficiently high tem- 

perature, and for power law potentials that satisfy the 
same condition as systematicity. We have 

Theorem. For power law potential traps (with power 
7 and dimension d), and for sufficiently high tempera- 
ture, the capacity of fermions is better than that of bosons 
when I/7+ 1/2 > 1/d. 

This includes, e.g., the harmonic trap in 2D and 3D, 
the 3D rectangular box, and the 3D spherical box. The 
theorem holds for quite moderate T, since we work up to 
the order (N/S\) 3 . Numerical simulations show indeed 
that it holds also for low temperatures, as seen, e.g., in 
Fig. [21 for 100 spinless fermions and same number of 




FIG. 2: We compare the channel capacity for the case of 100 
(noninteracting) spinless fermions (upper curve) with that of 
100 bosons (lower curve) . The trap is a 3D box with pbc. For 
the lower curve, the horizontal axis is T/T c , while for the up- 
per curve, it is T/T f , where T f = (h 2 /2mfcs)(67r 2 N/ gL 3 ) 2/3 
is the Fermi temperature. 



spinless bosons trapped in a 3D box with pbc. 

Note that we have found that the capacities for a 3D 
box without pbc are lower than those with pbc, both for 
bosons and fermions. The capacity for a harmonic trap, 
with the same characteristic length scale, has a higher ca- 
pacity (both for bosons as well as for fermions) than that 
for a 3D box with pbc. More importantly, the capacities 
for the case of fermions do not show any signatures of 
criticality, as expected. 

Interacting fermions. Until now, we have been deal- 
ing with the case of noninteracting bosons and fermions. 
Although a system of noninteracting massive bosons ex- 
hibits condensation, this is not the case for noninteract- 
ing fermions. A system of interacting fermions, however, 
can exhibit Cooper pairing, and, consequently, superfluid 
BCS transition (see e.g. 112]). It is therefore interesting 
to see whether such a "condensation" can be observed in 
the capacity of a channel transmitting trapped interact- 
ing fermions. We consider here a 3D box (of volume 
L 3 ) with pbc, within which iV fermions are trapped. 
The fermions behave like an ideal Fermi-Dirac gas, ex- 
cept when pairs of them with equal and opposite mo- 
mentum, and opposite spin components have kinetic en- 
ergy within an interval Ae on either side of the Fermi 
surface. In that case, the pairs experience a weak at- 
traction. The Hamiltonian can then be written as 7i = 

E fc -%(4, T a M + 4,i a -?J + ^^AA^-mv 

where t? = h 2 k 2 /2m, whereas V? vanishes except when 
\fjL - h 2 k 2 /2m\ < Ae and |/x - h 2 l 2 /2m\ < Ae, in 
which case, r = ~Vo < 0. Using mean field ap- 
proximation, the average occupation number in this case 
turns out to be = (1 - (eg/£ , g)tanh( y 8^/2))/2, 

where eg = t % - fx, E % = ^l+\& k -\ 2 . Here A £ = 
A(T), when |ed < Ae, and vanishing otherwise. This 
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FIG. 3: Channel capacity of 100 interacting fermions in a 
3D box trap with pbc, against T/Tf. The capacity is ini- 
tially convex, and then becomes concave, for higher T's, as 
illustrated in the inset. 



is the so-called "gap", given by the equation, Af = 
Vq J2j:(Aj:/Ej:) tanh(/3£^/2), where the summation runs 
only for an energy interval Ae about the Fermi surface. 
Note that the occupation numbers in this case are differ- 
ent from the case of ideal (noninteracting) fermions. Us- 
ing the gap equation and the constraint J2k — N on 
the total number of fermions in the trap, we can find the 
channel capacity by replacing n{ by iVg in C* E N . As we 
see in Fig. the channel capacity again changes its be- 
havior qualitatively, indicating the onset of the superfluid 
BCS transition. In Fig. [3] we plot the channel capacity 
against T/Tf, where Tf is the Fermi temperature for the 
case of noninteracting fermions. This is for convenience, 
as the thermodynamical transition temperature (for in- 
teracting fermions) has an exponentially decaying factor, 
which renders it inconvenient for our purposes. Also we 
chose Ae = and Vn = 10 -6 -Vj in the figure. 

In this Letter, we have used the GCE, which has ap- 
peared due to the dual constraints of average energy and 
average particle number. For a fixed number of particles, 
and retaining the average energy constraint, we are led 
to CE. In the thermodynamical limit, the average occu- 
pation numbers are the same for CE and GCE. For finite 
N , exact calculations for CE are difficult. However, dif- 
ferent approximate methods (see e.g. 0|) reveal that 
the average occupation numbers of CE are more uniform 
as compared to GCE, so that the former ensemble has 
larger capacity. However, the difference is marginal. 

The channels that we have considered in this Letter 
are noiseless. A simple, but physically important, model 
of noise is the Gaussian noise acting similarly on each 
mode, resulting in an effective increase of temperature in 
the channel. So for a given temperature, to accomodate 
the average energy constraint, we must start with a lower 
temperature than that in the noiseless case, leading to a 
decrease in capacity. The lower capacity in this particular 
noisy case can be read off from the figures of the noiseless 



one after finding the temperature difference. 

To conclude, we have considered the classical capacities 
of noiseless quantum channels carrying a finite average 
number of massive bosons or fermions. We have shown 
that the capacities are attained on the grand-canonical 
ensemble of statistical mechanics. Capacity of a channel 
carrying bosons indicates the onset of Bose-Einstein con- 
densation, by changing its behavior from being concave 
to convex with respect to the temperature, at the tran- 
sition point. Also the signature of the onset of Bardeen- 
Cooper-Schrieffer transition can be observed for weakly 
interacting fermions. We show analytically that for non- 
interacting particles, fermionic channels are better than 
the bosonic ones, in a wide variety of cases. 
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